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Multiply-connected Bose-Einstein condensed alkali gases: 
Current-carrying states and their decay 
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The ability to support metastable current-carrying states in multiply-connected settings is one 
of the prime signatures of superfluidity. Such states are investigated theoretically for the case of 
trapped Bose condensed alkali gases, particularly with regard to the rate at which they decay via 
thermal fluctuations. The lifetimes of metastable currents can be either longer or shorter than 
experimental time-scales. A scheme for the experimental detection of metastable states is sketched. 
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Multiply-connected superfluid and superconducting 
systems can support states in which a persistent macro- 
scopic particle-current flows. While not truly eternal, 
these states can have extraordinarily long life-times, their 
decay requiring the occurrence of certain relatively infre- 
quent but nevertheless topologically accessible (quantum 
or thermal) collective fluctuations (l[jj]. With the many 
considerable successes and rapid progress in the exper- 
imental exploration of Bose-Einstein condensed (BEG) 
alkali gas systems it seems reasonable to anticipate 
that multiply-connected settings for BEG will soon be- 
come available, thus allowing superfluid properties such 
as persistent currents to be sought. The purpose of the 
present Paper is to address, theoretically, the ability of 
BEG alkali gas systems in multiply-connected settings to 
support metastable current-carrying states, and to ad- 
dress the stability and decay of such states via thermal 
fluctuations. Gomplementary work by Rokhsar 1^ ad- 
dresses related questions regarding the creation of these 
states, and their stability. 

We adopt a phenomenological description in which we 
characterize the state of the BEG system by a macro- 
scopic wavefunction ^(r), in terms of which the conden- 
sate density n and current density j are given by 
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The free energy T of the state is given by the Gross- 
Pitaevskii form 
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where m is the mass of an individual atom, V^(r) is an 
effective external potential describing the magnetic and 
optical confinement of the atoms, and g (= Anh^a/m) 
represents the interatomic interaction, with a being an 
effective scattering length. For the sake of simplicity 
we neglect any possible effects due to spin. This de- 
scription is appropriate for analyzing the behavior of a 
BEG system at chemical potential fi and temperature 



T 0. In order that the condensate be able to undergo 
the free-energy (and angular-momentum) changing fluc- 
tuation necessary for current-dissipation the condensate 
must not be isolated. Therefore we restrict the system 
to be at temperatures not far below the critical tempera- 
ture Tc, in which case the non-condensed atoms serve to 
provide an energy and angular momentum reservoir. 




FIG. 1. Envisaged geometry of a trap supporting 
metastable current carrying BEG states. The condensate 
healing length ^ is regarded as being small, compared with 
the circumference of the torus L (= 2nf), but larger than its 
thickness R. 



As our aim is to address multiply-connected systems, 
we consider trap potentials V{r) that confine the gas to a 
cylindrically symmetric toroidal region (Fig. Hence, 
V depends only on r and z, where {r, (p, z} are the usual 
cylindrical polar coordinates. Moreover, we restrict our 
attention to systems in which the circumference of the 
torus L (= 27rf) is considerably greater than the conden- 
sate healing length ^ {fi^ /mgn)^^'^ , where n is related 
to the maximum particle-density], and the thickness of 
the torus R is comparable to or smaller than ^ This 
corresponds to a regime of low condensate-density; hence, 
the Thomas-Fermi approximation M is not applicable. 
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Traps having these gross features should be achievable 
by the use of magnetic and optical forces jl^] . There are 
two main reasons for considering this setting: (i) there 
would be no locally-stable current-carrying states if L 
were comparable to or smaller than ^; (ii) for thicker 
samples (i.e. R > (,), the relevant dissipative processes by 
which the persistent current decays become significantly 
more complicated (ultimately involving the nucleation of 
vortex rings) . 

We shall be concerned with events in which the system 
decays from some metastable current-carrying state 
(which is a local minimum of J-) to a lower-energy (and 
typically more stable) state via a thermal fluctuation . 
The current decays through a dissipative process during 
which the condensate density shrinks in magnitude over 
a region whose length is comparable to ^. Dynamically, 
one can envisage this process as occurring via the pas- 
sage of a vortex across the sample: a free-energy barrier 
must be overcome for this event to occur. The height of 
this barrier SF is given by the difference between the free 
energy of (the metastable state) and that of the tran- 
sition state ^Pt, i.e. the lowest possible free-energy high 
point en-route through configuration space between the 
initial and final metastable states. This thermally acti- 
vated process should occur at a rate Uoe"^-'^ , where, as 
we shall discuss later, the attempt frequency Wo does not 
contribute significantly to the temperature-dependence 
of the rate. 

In order to calculate the barrier heights, we first iden- 
tify the collection of metastable current-carrying states 
{^'m} and the relevant states {4't} for transitions be- 
tween them. Both families of states are stationary points 
of and therefore satisfy the time-independent Gross- 
Pitaevskii equation 

^ = -^^^^ + {V{v)-^i)^ + gm''^^Q. (3) 
subject to periodic boundary conditions in the coordinate 

<!>■ _ 

To address Eq. (|^) we introduce a complete orthonor- 
mal set of "transverse" eigenfunctions H^{r, z) and the 
associated energy eigenvalues A^, labeled by the "chan- 
nel" index v, which solve the eigenproblem 

- TT- (r-^drrdr + dl) + V{r, z) = KH,. (4) 
We then expand in terms of these eigenfunctions: 

^{r,^,z) = Y,FA^)HAr,z). (5) 

By inserting this expansion into Eq. (^ and making use 
of the orthogonality conditions for i?,^, we arrive at the 
following set of nonlinear coupled ordinary differential 
equations for the "longitudinal" wave functions F^i^c^): 

F: + a.F, -f3 F* F,^ F,^ = C., (6) 

Ul 1/2 1/3 



where primes denote derivatives with respect to 0, and 



the coefficients are defined by 

a^, = 2mf^ {ji — \y) /It?' , (7a) 

P = 2mf^g/h^, (7b) 

r;;^:;^ = [dzfrdrH, H,, H,, iJ,3 . (7c) 



The terms C^, given in footnote ||T^], are functions of the 
{Fjy} and are negligible when the transverse extent of the 
condensate R is small compared with the circumference 
of the torus L, scaling as [R/L)"^. 

The physical condition imposed above, viz., that ^ be 
much larger than the torus thickness i?, enforces the con- 
dition Ao < /i < Aj/^o- Hence, except iov v — we have 
< 0. In practice, is expected to be quite large, 
scaling as (i/i?)^. 

In the low-density limit the u — Q channel should domi- 
nate, all other channels being occupied weakly. We incor- 
porate this notion by introducing a book-keeping param- 
eter A into Eq. (|6|): for ^ we make the replacement 
/3 A/3. It can be verified that the non-leading terms in 
A may be neglected for the present purposes. (Incorpo- 
rating their effects is straightforward, if tedious.) Simi- 
larly, one can also incorporate the effects of the Ci, [F] . 
Hence, we find that the relevant states {^E'm} and {^t} 
can adequately be described in terms of Fo- 

Within the approximation scheme just outlined, the 
uniform current-carrying states have the form = 
fmc'^-'Ho, with 

fl ^ Nj27r - (a - n^)//3r, (8a) 
Sm = rimcl), (8b) 

for integral n^. For the sake of brevity, we now write 
a and T in place of ao and r°°. At the stated level 
of approximation, is the number of condensed par- 
ticles in the metastable state and F"^ is R^L, i.e., the 
volume occupied by the condensate. By considering the 
second variation of J- it can be readily shown that these 
states are local minima (and hence metastable) provided 
nC < I, where Cm = (47r//3F7Vm)^/^ ~ ^tt^/L is the di- 
mensionless coherence length. (This limit on the maxi- 
mum stable value of is the same as one would find 
using Landau's criterion for the critical velocity.) 
The transition states '^t — ft e''^' Ho are given by 

/2 = (iVt/2^) (1 - A2sech^(A(/./Ct)) , (9a) 
f^d^St^iNt/2Tr)nt. (9b) 

Far from a region of length ^, the amplitude ft is con- 
stant {ft ^ Nt/2-K) and the phase St winds uniformly 
{St ^ nt(j)). The coefficients in Eq. (|^) appear simplest 
when expressed in terms of the dimensionless coherence 
length Ct = (47r//?F7Vt)'/': 
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7Vt/27r ={a- n^^)/^3^, 

m ^ n~ TT^^ cos~-^(nt Ct), 
A2 = l-(ntCt)'- 



(10a) 
(10b) 
(10c) 



As N„j and iVt differ only by quantities of order ^/L, 
either of them may be used to characterize the num- 
ber of condensed particles. The transition states must 
have the property that they are saddle points of J- with 
only one direction of negative curvature. (This unsta- 
ble direction is the relevant reaction coordinate.) It can 
readily be shown that the states in Eq. (^ satisfy this 
condition as long as ^ > i?, or equivalently /i < X^^a- 
Thus, our approximation scheme for a BEG in a three- 
dimensional trap reduces the problem precisely to the 
one-dimensional problem addressed by Little |^ , Langer 
and Ambegaokar ||^ , and McCumber and Halperin Q . A 
useful by-product of the present approach is that it pro- 
vides a scheme for determining the intrinsic resistance 
of superconducting wires clad by normal-state materials 
(and thus having proximity-induced superconductivity). 

Having found the relevant states, we now calculate the 
free-energy barrier for dissipative fluctuations. It can be 
shown that states ^ satisfying Eq. (H) have a free-energy 
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Using this expression, along with Eq. (y) and Eq. (|9|), we 
find that 

5T = ]^5To [A (2 + (nt Ct)') - 3nt Ct cos-^n^ Ct)] , (12) 

where 5Tq is the long wavelength (i.e. rit 0) value of 
i.e., 
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We now develop order-of-magnitude estimates for the 
decay rates of metastable states via thermal fluctua- 
tions. Let us consider ^^Rb, for which the scattering 
length a is 5.8nm. We take a harmonic trapping po- 
tential V{r) = {l/2)mL0^[{r — f)^ + z^], whose ground 
state width y^Ti/moj can be identified with the width 
of the condensate R. To estimate we consider N 
non-interacting atoms in the potential V{r). By virtue 
of the geometry (i.e. R L) we can ignore the 
curvature of the torus, giving us a density of states 
p{E) = (4/3)(l/?itj)2(mLV27r2?i2)i/2£:3/2^ Integrat- 
ing this with the Bose occupation factor reveals that 
Tc « 1.28{h'^/m){N/R'^Lf/^. For example, if we assume 
that N « 10*^, Nt w 2.5 x 10^, R « l^im, and L « 100/im 
then we find 6fo /^b = 3.2fiK, and = 0.28^K. The 
barrier height is sensitive to changes in A't , and can there- 
fore be manipulated by heating or cooling the sample. 

The Arrhenius formula for the decay rate in terms of 
the barrier height is F « WoC"'''^^'^^. The attempt fre- 
quency ujo can be estimated by using the value of the 



microscopic relaxation time r, together with the assump- 
tion that each coherence volume in the sample fiuctuates 
independently [Q. A realistic estimate for t is the clas- 
sical collision time for a dilute gas [i.e. ~ anv ^ 
a'^{N/V){kBT/my^^ ^ 5 x lO'^ Hz], giving lifetimes for 
the metastable states on the order of seconds. Even be- 
yond the limits of validity of our calculation, one ex- 
pects SJ- to be a monotonically increasing function of 
the density. Hence, the barriers can be extremely large 
at low temperatures, allowing a continuous tuning of 
the metastable state lifetime from microseconds to times 
longer than the lifetime of the condensate. 

We now discuss two of the issues necessary for the ex- 
perimental testing of the predictions presented in this Pa- 
per. We have been considering the decay of metastable 
states, but have not yet addressed the issue of how to cre- 
ate them. Various approaches to creating a metastable 
current carrying state that rely on the superfluid proper- 
ties of the condensate have been discussed in detail by 
Rokhsar Q. Another technique, which does not rely 
on the superfluid nature of the condensate, takes ad- 
vantage of the spatial separation of the condensed and 
non-condensed atoms. One imagines starting the whole 
system rotating (for instance by using a rotating non- 
axisymmetric field) then applying a localized perturba- 
tion (such as from a sharply focused laser) that stops the 
thermal atoms but leaves the condensate rotating. 

The second experimental matter to be addressed is 
the detection of metastable current-carrying states. Per- 
haps the least difficult scheme would make use of present 
phonon imaging techniques [l4| . The experimental con- 
figuration could be as follows: A pulse of laser light 
generates a local rarefaction of the condensate, which 
then travels as two waves, one moving clockwise, the 
other counterclockwise. By non-destructive imaging 
techniques one might then observe where the two waves 
meet, which gives the velocity of the metastable super- 
current. This is only feasible if the speed of sound c 
is comparable to the velocity v with which the conden- 
sate moves around the annulus. Linearizing Eq. (|^) gives 
c = {gT / 271 my ^"^ w 1.2mm/s, which is only 30 times 
greater than v — h/mf w 46/im/s. 
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